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We analyze s-wave pion-nucleon scattering in a unitarized chiral effective Lagrangian including all dimension two contact terms. 
We find that both the 5n(1535) and the Sn(1650) are dynamically generated, but the S3i(1620) is not. We further discuss the 
structure of these dynamically generated resonances. 
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1. Introduction 

Pion-nucleon scattering has traditionally been the premier 
reaction to study the resonance excitations of the nucleon. In 
particular, in the S 1 1 partial wave, one finds two close-by res- 
onances at 1535 and 1650 MeV, which overlap within their 
widths of about 100 MeV. It was pointed out early in the frame- 
work of unitarized coupled-channel chiral perturbation theory 
fl] that this resonance might not be a three-quark (pre-existing) 
resonance but rather is generated by strong channel couplings, 
with a dominant KY. - KA component in its wave function. This 
analysis was extended in Ref. [2], where within certain approxi- 
mations the effects of 3-body nnN channels were also included. 
Further progress was madeinRef. [3], where the S n phase shift 
was fitted from threshold to about -\fs 2 GeV together with 
cross section data for n~p — > rjn and n~p — > K°A in the respec- 
tive threshold regions. This led to a satisfactory description of 
the S 1 1 phase and a reasonable description of the inelasticity 
up to the rjN threshold. Two poles were found corresponding 
to the 5 n(1535) and the 5 1 1 (1650) resonances together with 
a close-by unphysical pole on the first Riemann sheet. More 
recently, it was pointed out in a state-of-the-art unitary meson- 
exchange model that there is indeed strong resonance interfer- 
ence between the two Su resonances, as each of these reso- 
nances provides an energy-dependent background in the region 
of the other fl. 

In view of these developments and our attempts to construct 
a unitary and gauge-invariant model for Goldstone-boson pho- 
toproduction off nucleons based on coupled-channel unitarized 
chiral perturbation theory [5], we consider in this letter the two 
s-waves Su and Sn in pion-nucleon scattering. We work in 
the framework of a coupled-channel Bethe-Salpeter equation 
(BSE) including in the driving potential all local terms of sec- 
ond order in the chiral counting, thus going beyond the of- 
ten used approximation of simply including the leading order 



Weinberg-Tomozawa interaction. Further, we do not perform 
the often used on-shell approximation. Note that Kp scat- 
tering including such dimension two terms was already ana- 
lyzed in a framework equivalent to the on-shell approximation 
of the Bethe-Salpeter equation in Refs. J^rS. Our investiga- 
tion is restricted to center-of-mass energies below 1.8 GeV, as 
required for the future meson photoproduction studies. As we 
will show, both resonances in the S 1 1 partial wave are dynam- 
ically generated, even if the scattering data are fitted only up 
to V* — 1-56 GeV. Quite in contrast, the S 31 (1620) resonance 
is not generated by the coupled-channel dynamics. We also 
analyze the structure of the dynamically generated resonances 
as revealed through their coupling to the various meson-baryon 
channels. 



2. Formalism 

We consider the process of meson-baryon scattering at low 
energies. The s-wave interaction near the thresholds is domi- 
nated by the Weinberg-Tomozawa contact term, derived from 
the effective chiral Lagrangian 



4'i = (Bdy^DX - m )B) 
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where (. . .) denotes the trace in flavor space, D^B := d^B + 
\[W , d^u], B], mo is the baryon octet mass in the chiral SU(3) 
limit, and D, F are the axial coupling constants. The relevant 
degrees of freedom are the Goldstone bosons described by the 
traceless meson matrix U, 
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where Fq is the meson decay constant in the chiral limit, and 
the low-lying baryons are collected in a traceless matrix 



B = 



V2 



A 
V6 






P * 






A 


n 






V6 




no 




vs v 



(3) 



We set external currents to zero except for the scalar one, 
which is set equal to the quark mass matrix, s — M : — 
diag{m u , mj,m s ). We furthermore use 



m 2 := U , uf := iu'd^u - iudPu* , 
^± := m'^m 1 " + ux T u , ^ := 2B s , 



where the constant Bq is related to the quark condensate in the 
chiral limit. 

The Weinberg-Tomozawa contact term mentioned above 
stems from the covariant derivative D^B , and is of first or- 
der in the chiral power counting. Most chiral unitary ap- 
proaches restrict their meson-baryon potential to this interac- 
tion, which generates the leading contribution to the s-wave 
scattering lengths. This approach has been remarkably success- 
ful in many cases, see e. g. Illlsil-Hsll. However, at first chiral or- 
der, there are also the Born graphs, describing the s-channel and 
M-channel exchanges of an intermediate nucleon. The full inclu- 
sion of these graphs in the driving term of the Bethe-Salpeter 
equation leads to conceptional and practical difficulties, which 
have not yet been solved to the best of our knowledge: (i) Iter- 
ation of the s-channel exchange Born graphs will generate var- 
ious contributions leading to a renormalization of the various 
baryon masses (and wave function renormalizations), which are 
usually set to their physical values in the loop functions of the 
chiral unitary approach. These contributions would thus have 
to be dropped. In view of a later application to photoproduc- 
tion, such a non-perturbative treatment of s-channel exchanges 
leads to complications with gauge invariance because the self- 
energies are linked (via a Ward-Takahashi identity) to the elec- 
tromagnetic baryon form factors, which would also have to be 
treated in a corresponding (non-perturbative) fashion, (ii) Iter- 
ation of the M-channel diagram, on the other hand, leads to all 
kinds of genuine multi-loop topologies, as there is no factor- 
ization into simple one-loop terms any more. The correspond- 
ing integral equation could only be solved numerically, e.g. by 
employing a Wick rotation and a four-momentum cutoff. Prob- 
lems with gauge invariance would also occur here. In the liter- 
ature, the M-channel Born diagrams were usually treated within 
some approximation which effectively reduced the solution of 
the BSE to products of one-loop terms, or included perturba- 
tively to guarantee a matching to ChPT amplitudes up to a given 
order, see e.g. [14]. All these approximations, however, destroy 
the exact correspondence of the individual terms in the solu- 
tion of the BSE to dimensionally regularized Feynman graphs, 
which is crucial in our approach to photoproduction. Therefore, 
we will approximate our interaction kernel by a sum of contact 
terms. To go beyond the simple Weinberg-Tomozawa potential, 
we shall include the full set of meson-baryon vertices from the 



second order chiral Lagrangian. These terms may lead to size- 
able corrections to the leading-order results, see e. g. the calcu- 
lation of NNLO corrections on meson-baryon scattering lengths 
within SU(3) ChPT 11511 . The pertinent Lagrangian density was 
first constructed in [ 16] and reads in its minimal form IU7I1 

4 2 i = b D , F (B\x + , B]J + b (BB)(x + ) 
+ b X i 2 (B[u^ [«", fl]J> + h(B[ Ull , [if, B}\) + b 4 (BB)(uX) 
+ Mb5/6<fic^ v [[ MjU , K v ], B\) + ib-iiB^u^iuyB) 

ibva/.— 



(4) + j^((Bf [u„, {My, [D\ B]}}) + (Bf[D v , [u\ {«„ B}}])) 
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(5) 



with the bi the pertinent dimension-two low energy constants 
(LECs). The LECs £>o,d,f are the so-called symmetry breakers 
while the bid - 1, . . . , 1 1) are referred to as dynamical LECs. 

The strict perturbative chiral expansion is only applicable at 
low energies. Moreover, it certainly fails in the vicinity of res- 
onances. The purpose of the present work is the extension of 
the range of applicability of the low-energy effective theory by 
means of a coupled channel Bethe-Salpeter equation (BSE). In- 
troduced in 111 (811 it has been proven to be very useful both in the 
purely mesonic and in the meson-baryon sector CI HUB. In 
contrast to perturbative calculations this approach implements 
two-body unitarity exactly and in principle allows to generate 
resonances dynamically. Due to the exact correspondence of 
the Bethe-Salpeter scattering amplitude with an infinite sum of 
dimensionally regularized Feynman graphs, we can use our so- 
lution of the BSE as an extended vertex in a model amplitude 
for meson photoproduction and arrive at a natural and straight- 
forward way to implement gauge invariance in a chiral unitary 
framework (for details on the construction principles, see J5|). 

In this section we collect the necessary formalism of the 
Bethe-Salpeter approach. We denote the in- and outgoing me- 
son momenta by q\ and respectively. Moreover the overall 
four-momentum is given by p — q\ + p\ — q 2 + pi, where p\ 
and p 2 are the momenta of in- and out-going baryon, respec- 
tively. For the meson-baryon scattering amplitude T^,^; p) 
and chiral potential V^,^;/?) the integral equation to solve 
reads 

T(<t 2 4vP) = VQt 2 ,ti,p)+ 

if ^V(4 2 ,kp)S(p-l)A(t)T0,<tvP), (6) 

where S and A represent the baryon (of mass m) and the me- 
son (of mass M) propagator, respectively, and are given by 
iS(f) = - m + ie) and iA(k) = i/(k 2 - M 2 + ie). The BSE 
is depicted in Fig.Q] 

So far we have suppressed the channel indices in the above 
formulas, however since we are dealing with coupled channels, 
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Figure 1 : Symbolical representation of the Bethe-Salpeter equation. Here the 
square and the circle represent the potential V and the scattering amplitude T, 
respectively. 



T, V, S and A are matrices in channel space (the propagators are 
represented by diagonal matrices). In view of a later application 
to photoproduction off protons, we restrict ourselves to meson- 
baryon channels with strangeness S = and electric charge 
Q = +1 . This leaves us with the following channels: 
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pn 



, AK + , Z°K + , Z + K° . 



(7) 



Now let us specify the interaction kernel to be iterated by means 
of Eq. (|6]l. As explained above, we only include the contact- 
term contributions from Xig and jjQ and omit the Born terms. 
To our knowledge this is the first time these NLO corrections of 
the chiral potential are included and unitarized within the full 
relativistic BSE, without making use of the on-shell approxima- 
tion or s-wave projection of the chiral potential, so that also a 
p-wave is iterated. Separating the momentum space from chan- 
nel space structures the chiral potential considered here takes 
the form: 

+ Ai A {qi-q 2 )+A il [qli,qf 1 \+A M {qi-q 1 ) 
+ A & i , (qk(qi ■ p) + ql\ (q 2 ■ p)), (8) 

where the first matrix only depends on the meson decay con- 
stants F„, Fk, F„, whereas A14, A^i, A%\\ and Am also contain 



the NLO LECs as specified in Appendix A In going from the 
Lagrangian (0 to the above vertex rule, we have left out some 
terms which are formally of third chiral order. 

The loop diagrams appearing in the BSE Eq. (0 are in gen- 
eral divergent and require renormalization. In case of a strict 
chiral perturbation expansion, the terms can be renormalized 
in a quite straightforward way, order by order, including at a 
given order of the calculation all the counterterms absorbing 
the loop divergencies. On the other hand the treatment of the 
divergencies of the BSE is known to be a complicated issue, 
see e.g. dUGji. Although the unitarization of the chiral po- 
tential provides us with large benefits regarding dynamically 
generated resonances, it relies on approximations of the kernel, 
which destroy some fundamental features of quantum field the- 
ory, such as crossing symmetry. 

There are various ways to treat the divergent integrals and the 
large baryon mass scale appearing. Without going into details 
here, we preserve the analytic structure of the loop integrals by 
utilizing dimensional regularization and just replacing the di- 
vergent part by a subtraction constant. The purely baryonic in- 
tegrals are set to zero from the beginning. Thus, our treatment 
of the loop integrals is, in effect, similar to the EOMS regu- 
larization scheme advocated in I20I1 . As it was argued in Q5|] 



it is not possible to express the terms necessary to absorb the 
divergencies in the BSE as counterterms derived from a local 
Lagrangian. However it is possible to alter the loop integrals 
in the solution of the BSE in a way that is in principle equiv- 
alent to a proper modification of the chiral potential itself (for 
an explicit demonstration, see App. F of 12111 '). In this spirit 
we apply the usual MS subtraction scheme, keeping in mind 
that the modified loop integrals are still scale-dependent. This 
regularization scale (p) dependence would be canceled by the 
corresponding scale dependence of higher-order contact terms 
in the perturbative approach, but in our nonperturbative frame- 
work, the scale /j is used as a fitting parameter, reflecting the 
influence of higher order terms not included in our potential. 
Note that in [l9|, the 12 loop integrals (4 for each meson- 
baryon, meson and baryon case) appearing there, gave rise to 
12 finite subtraction constants, which were then also used as 
fitting parameters of their approach. 

Having specified the kernel we are now ready to solve the 
Bethe-Salpeter equation. Given the structure of the kernel, its 
iteration via the BSE induces the following form of the scatter- 
ing amplitude, 
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(9) 



1=1 



where the coefficients T, are 6 x 6 matrices in channel space, 
which only depend on the center-of-mass energy after fixing 
the LECs, and N := (qli, pq/\, qlifqix, qkql\, pq!\{qi - p),qi\(qi- p), 

qfi(q\ ■ p), qkq/i, (qi ■ p)(qi ■ p), p(q\ ■ p)(qi ■ p), (q\ ■ p), 
f(qi ■ p), (qi ■ q\), piqi • <?i ), qkp, qk, p(qi ■ p), (qi ■ p), l, p) is a 

vector in the 20-dimensional space of invariant structures. Note 
that the scalar products are listed here as independent structures 
because we include the full off-shell dependence of the chiral 
potential in the BSE, which prevents us from writing them as 
simple functions of the Mandelstam variables s and t. 

On the other hand the above decomposition allows us to pull 
the coefficients T, out of the loop-integral in Eq. Then 
these are fully determined by the solution of a linear system of 
equations in the space of invariant structures: 
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(U = i, 



,20). 



(10) 



where the V, are coefficients of the chiral potential with respect 
to the invariant structures defined above and X is a 20 x 20 
matrix. The latter connects different structures of the space of 
invariant structures via loop integrations on the r.h.s. of Eq. (0. 
Once the BSE has been solved, we can of course set the external 
four-momenta on their mass shells, leaving us with only two 
independent structures for the on-shell amplitude, i.e. 1 and f. 

3. Results and discussion 

Throughout the present work we use the following numerical 
values (in GeV) for the masses and the meson decay constants: 
F„ = F,,/l.3 = 0.0924, F K = 0.113, M„ a = 0.135, M K+ = 
0.1396, Mr, = 0.5478, M K+ = 0.4937, M K o = 0.4977, m p = 
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Figure 2: Real and imaginary part of the S\ \ partial wave amplitude compared 
with the SAID-data (WI08-analysis). Full curves correspond to the best fit, the 
dashed ones to fits with slightly worse x\r The bold vertical line limits the 
region of the fit, where in the non-fit region single energy values are taken from 
the SAID-data. 

0.9383, m n = 0.9396, m A = 1.1157, w E o = 1.1926 and m z + = 
1.1894. The baryon mass in the chiral limit, mo in Eq. (0, can 
be fixed to 1 GeV without loss of generality, as any other value 
only amounts to a rescaling of the unknown LECs. 

There are 17 free parameters in the present approach, given 
by the 14 LECs, as well as three subtraction constants for the 
regularized loop integrals, corresponding to the logarithms of 
the undetermined regularization scales (in GeV), i.e. log(jj.„), 
log(/i^) and log(yu^). Here we take the regularization scale of 
each channel to be fixed by the respective meson, i.e. in addi- 
tion to n„ N =: fi„ and fx nN =: fi n , we take n K1 _ = jj. KA =: fi K . The 
latter constraint appears to be natural in view of our forthcom- 
ing work on meson photoproduction, where loops are present 
in which a photon-induced A — > S° transition occurs. 

For the fits, we consider experimental data for s-wave nN 
scattering up to W = 1 .56 GeV, i.e. partial wave amplitudes 5 n 
and 531 (both real and imaginary parts) provided by the SAID- 
program at GWU, see |22ll. Comparing an earlier analysis by 
the Karlsruhe group 02311 to the current one, we assign for the 
energies below W = 1.28 GeV an absolute systematic error 
of 0.005 and for higher energies an error of 0.030 to the partial 
wave amplitudes. To some extent this is in agreement with error 
estimates done in which are motivated by the expectation 
of pronounced three-body effects above the kkN threshold. For 
the best fit, found using the MINUIT library, with a^ of = 1 .23 
we obtain the following parameter set (all in GeV -1 ) 



log(^) = 


+0.924, 


ft 4 = 


-0.215, 


bw = 


+ 1.920, 


logoff) = 


+0.581, 


b 5 = 


-0.963, 


b u = 


-0.919, 


log(yU^) = 


-0.218, 


b 6 = 


+0.218, 


b = 


-0.768, 


ft, = 


-0.082, 


bi = 


-1.266, 


b D = 


+0.641, 


ft 2 = 


-0.118, 


H = 


+0.609, 


b F = 


-0.098, 


^3 = 


-1.890, 


b 9 = 


-0.633. 








1200 1400 1600 1800 2000 

Wcms [MeV] 

Figure 3: Real and imaginary part of the S31 partial wave amplitude compared 
with the SAID-data (WI08-analysis). Full curves correspond to the best fit, the 
dashed ones to fits with slightly worse x\ f The bold vertical line limits the 
region of the fit, where in the non-fit region single energy values are taken from 
the SAID-data. 

All parameters are of natural size and LECs agree with the esti- 
mates from the SU(3) to SU(2) matching relations provided in 
lfl5ll . However we are only able to estimate the computational 
errors on the above parameters within the MIGRAD (MINUIT) 
minimization procedure, which appear to be negligible. 

In Figs.|2]and[3]we present the result of our approach for the 
5 11 and 531 partial waves. As already seen in earlier publica- 
tions on the BSE approach with leading order chiral potential 
fl, the low-energy region (e.g. yfs < 1.4 GeV) is reproduced 
for both isospin 3/2 and 1/2 reasonably well. For the two s- 
wave scattering lengths, we obtain 01/2 = 145.8 x 10 _3 /M ff + 
and 03/2 = -91.6 x 10 _3 /M ff +, to be compared with the direct 
extraction of these scattering lengths from the GWU solution, 
a V2 = (174.7 + 2.2) x lO" 3 /^ and a 3 /2 = (-89.4 ± 1.7) x 
lO^/Myr+Q The theoretically cleanest determination of these 
observables stems from the analysis of pionic hydrogen and pi- 
onic deuterium data based on effective field theory M24ll . fli/2 = 
(179.9 ± 3.6) x lO^/M^ and a y2 = (-78.5 ± 3.2) X lO^/M^. 
The description of the nN amplitude at low energies will cer- 
tainly be improved by a more complete treatment of the Born 
terms, which is beyond the scope of this Letter. One might also 
think about constraining the well-known pion-nucleon scatter- 
ing lengths, e.g. by adopting a matching procedure to the per- 
turbative expansion. However, since we did not put a special 
weight on the threshold region in our fits, and the overall de- 
scription of the partial waves seems to work well over a rather 
broad energy range, we regard the obtained results for the scat- 
tering lengths as satisfactory. 

Moreover, and more importantly, within the fit region we re- 
produce the 5n(1535), without any use of explicit vector me- 
son resonances or even taking into account the nnN channels as 
for example in [2]. At the same time the 5 31 (1620) resonance is 



'We thank Ron Workman for providing us with these values. 
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Figure 4: (222-1 1 1) Riemann sheet of the .v-plane. The five-star and the six-star 
correspond to the values obtained in Ref. |4] and Ref. [19], respectively, dots 
represent results of phenomenological models listed in I25fl . 




2.5 2.55 2.6 2.65 2.7 2.75 2.8 
Re s 



Figure 5: (2222-1 1) Riemann sheet of the s-plane. The five-star and the six-star 
correspond to the values obtained in Ref. [4] and Ref. [19], respectively, dots 
represent results of phenomenological models listed in |2al . 



not reproduced by our approach, which is in agreement with the 
current state of knowledge that the first S31 resonance does not 
have a prominent dynamically generated component. To em- 
phasize this we exclude the data on S 31 and recalculate the^ of 
for the above parameter set, we end up with^ of (5 n) = 0.59. 

At this point one realizes an even more interesting fact: 
After fixing the S u partial wave in the energy region up to 
sfs — 1 .560 Ge V every curve with minimized )£. possesses 
a second structure between KA and KI. threshold. Obviously 
this corresponds to the well-known S n (1650) resonance and 
is predicted here only by demanding a good description in the 
low-energy and the first resonance region. To some extent this 
is in agreement with Ref. |3J|, where the S n(1650) was repro- 
duced in the fit of the phase shifts and inelasticities for the full 
region of 1.077 < yfs/GeV < 1.946. While only the leading 
order chiral potential was considered there, the authors intro- 
duced additional parameters appearing for every loop integral. 
Apparently these parameters contain some of the information 
that has to be attributed to neglected terms of higher order in 
the chiral potential. Additionally, in contrast to our approach 
this method does not allow to identify the higher partial waves 
than the s-wave, which might become important for higher en- 
ergies as emphasized in |5|]. 

In Figs.|4]and|5]we present the modulus of the analytic con- 
tinuation of T}_h into the complex s-plane. In Fig. [4] two poles 
appear on the (222-111) Riemann sheet, which labels the un- 
physical Riemann sheet connected to the physical (scattering) 
axis in the energy region between the third and fourth threshold, 
i.e. (M n + mjv) 2 < s < {Mk + m^) 2 . For the position of the two 



poles we extract: 

W1535 = (1.506- 0.140 OGeV, 

W I650 = (1-692 - 0.046 GeV. (12) 

Choosing the (2222-11) Riemann sheet, i.e. the unphysical 
sheet reached by analytic continuation from the region (Mk + 
ika) 2 < s < {Mk + '«l) 2 , see Fig.|5j we obtain one single pole 
structure, which is located at 

W I650 = (1-682 - 0.042 i) GeV. (13) 

We conclude that the S n(1650) can also be described as a dy- 
namically generated resonance, just like the S n(1535). 

Clearly the uncertainty of our predictions grows with increas- 
ing energy. As a consequence of the sizeably increased com- 
puting time, when fitting the full amplitudes rather than the on- 
shell approximations to them, we are not able to perform a full 
error analysis as e.g. done in Ref. [8] for K p scattering. Still, 
we are able to get an indication of the error bands on the par- 
tial wave amplitudes. For this we present the second, third and 
fourth best fits in Figs. [2] and [3] as dashed lines. However the 
error analysis deserves further studies. 

It is further interesting to analyze the structure of these states. 
To do that, we consider the on-shell scattering matrix in the 
vicinity of the two poles, where it takes the form 

- — - , (14) 

1 S-S R 

with gi {gj) the complex coupling constant for the initial (final) 
transition of the meson-baryon system. For the S n(1535), we 
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obtain the following ordering 

\gAK^ > \g P / > ISX+K0| 2 - \g^\ 2 > \gz°K+\ 2 - Igpiflt d5) 

We remark that the inequalities between couplings to different 
nN and KY, channels are mostly due to Clebsch-Gordan co- 
efficients in the associated isospin decompositions. However, 
isospin symmetry is not exact in the present approach. We find 
that the largest component is the KK one and that the coupling 
to tjN is significantly bigger than the nN ones, in agreement 
with the empirical fact that the S n(1535) couples dominantly 
to TjN. The pattern for the S n(1650) looks different, 

\gl.*K»\ 2 > \gp/ > l^ + | 2 * \g m+ \ 2 > \g pn o\ 2 » \g AK +\ 2 , (16) 

i.e. for this resonance the KI. component is dominant and the 
KA one is completely negligible, which for instance is indi- 
cated by the fact that the pole associated with the S n(1650) is 
accompanied by a second one on a neighboring sheet, with al- 
most the same coordinates. As for the lower-lying resonance, 
the coupling to Nq is bigger than the one to Nn. 

4. Summary and outlook 

In this Letter, we have analyzed s-wave pion-nucleon scat- 
tering in coupled-channel unitarized chiral perturbation theory. 
The driving kernel includes all local interactions terms of first 
and second order from the chiral effective Lagrangian. We con- 
sider all two-body channels with strangeness zero and charge 
plus one, but do not include inelasticities generated from three- 
body Nnn states. The Bethe-Salpeter equation has been solved 
including the full off-shell dependence of the chiral potential. 
The parameters are fitted to the real and imaginary part of the 
S a and the S 31 partial waves for cms energy below 1.56 GeV. 
We show that both the Sn(1535) and the Sn(1650) are gen- 
erated dynamically, even though the fit range does only include 
the first resonance. We have also analyzed the structure of these 
states, which exhibit some marked differences as indicated by 
the couplings given in Eqs. ( 1151161 1. Quite differently, no reso- 
nance is generated in the S31 partial wave. We consider this an 
important step in our program of describing kaon photoproduc- 
tion from coupled-channel unitarized chiral perturbation theory. 
Clearly, in the future more work is needed to properly include 
the Born terms and to perform a systematic error analysis. 
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Appendix A. Couplings 

For the channel indices {b, j; i, a} corresponding to the pro- 
cess <piB a — > 0jBh the relevant coupling matrices read 

+ b 2 ((A M {A'\ [A', A a ]}) + (A b \X, [A j \ A a ]})) 

+ /><(</ \A' : . {A 1 , A"}}) + {A 1 " {A, {A ji , A a }})) + 2b 4 {A bi A a ){A n A)), 

A b f a = b 5 (A b V,A i ],A a ]) + b 6 (A"H[A^,A i ],A a }) 

+ bj^A^XA'A ) + (A b ^A i ){A a A Jf ))), 

+ b 9 ({A bi [A j \ {A\ A"}]) + {A b \A\ [A j \ A a }])) 
+ b l0 ((A h HA jf , {A 1 , A a }}} + </V, [A fl , A a }}}) 
+ 2b n {A b ^ A a ){A* A 1 )), 

A b f'° = --l-(b D ({A b \{A 3 \{M,A%A a }) 

+ {A b \[A t ,{M,A*)},A a }j) 
+ b F {(A b ^[{A n , {M, A 1 }}, A"]} + {A b \{A\ {At, A*}}, A a ])) 
+ 2b {{A h " A a ){[A* A^M))), 

where A denote the 3x3 channel matrices (e.g. <p = <p'A' for 
the physical meson fields), the F ,■ are the decay constants of 
the meson in the respective channel, and (...) denotes the trace 
in flavor space. Moreover, M. is obtained from the quark mass 
matrix M via the Gell-Mann Oakes Renner relations, and given 
in terms of the meson masses as follows, M = |diag(M 2 f+ - 
M| + M 2 ,,, M2„ - M 2 + + M 2 ,,, M 2 k+ + M\ a - M 2 ^) . 
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